Nuclear spin-lattice relaxation in ferrimagnetic clusters and chains: 
A contrast between zero and one dimensions 
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Motivated by ferrimagnetic oligonuclear and chain compounds synthesized by Caneschi et at, 
both of which consist of alternating manganese(II) ions and nitronyl-nitroxide radicals, we calculate 
the nuclear spin-lattice relaxation rate 1/Ti employing a recently developed modified spin- wave 
theory. 1/Ti as a function of temperature drastically varies with the location of probe nuclei in both 
clusters and chains, though the relaxation time scale is much larger in zero dimension than in one 
dimension. 1/Ti as a function of an applied field in long chains forms a striking contrast to that 
in finite clusters, diverging with decreasing field like inverse square root at low temperatures and 
logarithmically at high temperatures. 

PACS numbers: 75.50.Xx, 75.50.Gg, 76.50.+g, 75.30.Cr 



I. INTRODUCTION 

Clusters of metal ions [1] serve to test the validity 
of quantum mechanical approaches at the nanome- 
ter scale. Among others a dodecanuclear manganese 
complex of formula [Mni 2 Oi2(CH 3 COO)i 6 (H 2 0)4] 
[2] and an octanuclear iron complex of formula 
[Fe 8 (N 3 C 6 Hi5) 6 02(OH)i 2 ] 8 + [3], which are hereafter 
abbreviated as Mni 2 and Fes, respectively, have been 
attracting considerable interest in this context. In these 
clusters, magnetic hysteresis curves are independent of 
temperature and consist of equally separated steps at 
sufficiently low temperatures [4-7], suggesting quantum 
tunneling of the magnetization. Both clusters have a 
magnetic ground state of total spin 10 and are often 
treated as a rigid spin-10 object with an Ising-type 
magneto-crystalline anisotropy. Such a phenomenologi- 
cal interpretation is indeed successful for the magnetic 
relaxation at sufficiently low temperatures and weak 
fields, but it masks individual internal structures of var- 
ious magnetic clusters. Without a microscopic Hamil- 
tonian, we could not essentially distinguish the Mni2 
and Fes clusters. We cannot compare nanomagnets with 
bulk magnets in the same microscopic language until we 
describe them in terms of their constituent ion spins. 

Thus motivated, several authors made an attempt to 



estimate exchange interactions in the Mni2 and Fes clus- 
ters calculating the low-lying energy spectra [8-10], mag- 
netic susceptibilities [11], magnetization curves [12,13], 
inelastic neutron-scattering spectra [14], tunneling split- 
tings [15], and nuclear spin-lattice relaxation rates [16]. 
The 8-spin modeling of the Mni2 cluster [14,15] revealed 
the decisive role of multispin effects in magnetic tunnel- 
ing, while the extended ab initio calculation [10] demon- 
strated the relevance of metal-ligand orbital bybridiza- 
tion to intramolecular exchange interactions. However, 
the predicted magnetic structures are rather controver- 
sial for both clusters. The Mn 12 cluster consists of eight 
Mn 3+ ions of spin 2 and four Mn 4+ ions of spin |, while 
the Fe§ cluster contains eight Fe 3+ ions, both of which 
are coupled to each other through four kinds of exchange 
interactions (Fig. 1). Antisymmetric exchange interac- 
tions of the Dzyaloshinsky-Moriya type are also assumed 
to be relevant to the Mni2 cluster [14]. There is no estab- 
lished way of assigning the magnetic anisotropy to each 
ion site. Thus, complicated intracluster magnetic struc- 
tures block our microscopic understanding of mesoscopic 
magnetism. Nanomagnets such as Mni2 and Fes have in- 
deed been providing fascinating observations but are not 
necessarily suitable for a comparative study on zero- and 
one (or higher)-dimensional quantum magnetism. 
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FIG. 1. Schematic plot of Mni2 (a), Fes (b), and (MnNIT)6 (c). Symmetry-inequivalent sites are distinguishably drawn. 
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In order to reveal how paramagnetic spins grow 
into bulk magnets, we take a great interest in fer- 
rimagnctic ring clusters [17] and chains [18] syn- 
thesized by Caneschi et al, both of which consist 
of manganese hexafluoroacetylacetonates (hfac) and 
nitronyl nitroxidc radicals 2-i?-4,4,5,5-tetramethyl-4,5- 
dihydro-l_ff-imidazolyl-l-oxyl 3-oxide (NIT-i?) with R = 
phenyl (Ph), isopropyl (i-Pr), ethyl (Et), methyl (Me). 
The [Mn(hfac)NIT-Ph]6 cluster, hereafter abbreviated as 
(MnNIT)6, has antiferromagnetic exchange coupling be- 
tween the Mn 2+ ions of spin | and the radicals of spin 
i so as to exhibit a magnetic ground state of total spin 
12 (Fig. 1). The [Mn(hfac)NIT-i?] 00 chain, hereafter 
abbreviated as (MnNIT)^, may be regarded as a one- 
dimensional analog of (MnNIT)g. Their simple magnetic 
structures, describable within isotropic exchange Hamil- 
tonians [17,18], are suitable enough to compare oligonu- 
clear ferrimagnets with those of one dimension in their 
intrinsic features. 

The theoretical tool we employ here is a recently de- 
veloped modified spin-wave theory, which is quite use- 
ful in understanding thermal [19-21] as well as ground- 
state [22,23] properties of various one-dimensional fer- 
rimagnets. We inquire further into zero dimension 
and dynamic properties. Besides direct observations 
of resonant magnetization tunneling, nuclcar-magnetic- 
resonance (NMR) measurements have extensively been 
performed for cluster magnets. 1 H, 2 D, 13 C, and 55 Mn 
NMR in Mnis [24-28], : H and 2 D NMR in Fc 8 [29,30], 
and 1 H NMR in a hexanuclear copper complex [31] 
significantly contributed toward revealing the quantum 
dynamics of cluster magnets. As for the systems of 
Mn(hfac)NIT-i?, we may consider NMR measurements 
using as probes 1 H, 13 C, 19 F, and 55 Mn nuclei. We cal- 
culate the nuclear spin-lattice relaxation rate 1/Ti as a 
function of temperature, an applied field, and the loca- 
tion of probe nuclei. 



II. MODEL HAMILTONIAN 



of 7io to Hi, we obtain the diagonal spin-wave Hamilto- 
nian as 

k 

(2.3) 

with E s = -2SsJN + J2 i= i,o E i and w± ( fc ) = 
J2 i=1 iof (k) , where Ei and wf(k) are the 0(S l ) cor- 
rections to the ground-state energy and the dispersion 
relations, respectively. The dispersions of the linear spin 
waves, tof(k), and the corrections due to the interactions 
between them, cj^(fc), are, respectively, given by 

uf(k) = [uj k ± (S- s)]J=F5MBff , 



cu±(k) = -2(S + s)jr 1 



sin 2 (/c/2) 



jr 2 



[uj k ±(S- s)] , 



(2.4) 



(2.5) 



where u k = [(S - s) 2 + 4Ss sin 2 (k/2)} 1 / 2 , 
A = -(l/2A0E fe [l - (S + s)/uj k ], and T 2 = 
(l/iV)E fc (V5^M)cos 2 (fc/2). A is nothing but 
the quantum spin reduction (^/N)^2 i {ajai)T=o = 
(1/-/V) Ej(&|&i)T=o an d can analytically be evaluated as 
l/v / 31 + 1/v^T- 7/24 at N = 6, which is slightly larger 
than the N — * oo numerical estimate 0.106139 and thus 
suggests growing quantum fluctuations with decreasing 
system size. In Fig. 2 we plot the spin-wave dispersions 
together with quantum Monte Carlo calculations. The 
lower branch, which reduces the ground-state magneti- 
zation, is of ferromagnetic aspect exhibiting a quadratic 
dispersion at small momenta, whereas the upper branch, 
which enhances the ground-state magnetization, is of 
antiferromagnetic aspect being gapped from the ground 
state. The antiferromagnetic mode is remarkable for its 
0(S°) quantum correction. 



(MnNIT) 6 and (MnNIT)^ are both described by an 
isotropic spin-(S', s) Heisenberg Hamiltonian [17,18] 



n = j 



N 

i=i 



N 



Si ■ S i+ i) - gusH^S? + sf) , 
i=i 

(2.1) 



where S = | 



and we have set their g factors both 
equal to g. Introducing bosonic operators for the spin 
deviation in each sublattice via S? = (2S — a|a i ) 1 / 2 a i , 



Sf = S - aUi, s 



&t(2s-&t6 i )V2 ) 8 z 



b\bi, we 



expand the Hamiltonian with respect to l/S as 

H = -2SsJN + Hi+H Q + 0(S- 1 ) 1 (2.2) 

where we assume that O(S) — O(s) and Hi is the 0(S l ) 
contribution. Considering the perturbational treatment 
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FIG. 2. Dispersion relations of the linear (LSW) and in- 
teracting (ISW) spin waves, tj^(k) and tof(k) + Uq (k), at 
N = 6 and N = oo under zero field. Corresponding quantum 
Monte Carlo calculations (QMC) at N = 32 are also shown 
for reference. 
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FIG. 3. Magnetic-susceptibility mea- 

surements on [Mn(hfac)NIT-Ph] 6 at H = 0.5 T [15] and on 
[Mn(hfac)NIT-i-Pr]oo at H — [16] compared with quantum 
Monte Carlo calculations (QMC). They are in good agree- 
ment assuming J/k B ~ 370 K for [Mn(hfac)NIT-Ph] 6 and 
J/k B ~ 360 K for [Mn(hfac)NIT-i-Pr]oc. Linear (LMSW) 
and interacting (IMSW) modified spin-wave calculations at 
N = 00 are also shown. 



III. MODIFIED SPIN- WAVE FORMALISM 

In order to calculate thermodynamics, we consider 
modifying the spin- wave scheme [32], that is, control- 
ling the number of bosons induced thermally. For 
isotropic ferrimagnets, we claim that the thermal fluctu- 
ation should cancel the staggered magnetization as [19] 



k cr=± 



Uk 



(3.1) 



where n£ 



J2n±=o n<Tp k( n ~i n+ ) witn Pk(n~,n + ) be- 
ing the probability of n~ ferromagnetic and n + antifcrro- 
magnetic spin waves appearing in the fc-momentum state. 
Minimizing the free energy 



k <J=± 



n , n" 



(3.2) 



k „± 



with respect to Pk at each k under the condition (3.1) to- 
gether with the trivial constraints E n ± Pk(n~,n + ) = 1, 
we obtain the optimum distribution functions as 

(3.3) 



k eacp{[w ± (fc)-/i(S , + a)/wfc]/*BT}-l' 

with a Lagrange multiplier \x due to Eq. (3.1), where 
the dispersions w ± (fc) are restricted to u>f(k) or set to 
cof (k) + ui^ (k) according as we consider linear or inter- 
acting modified spin waves. In the present formalism, 
the dispersion relations are determined without modi- 
fying the original Hamiltonian (2.2) and then the La- 
grange multiplier is introduced so as to construct a reli- 
able thermodynamics, which is essentially different from 
the Takahashi scheme [32] and enables us to investigate 



much wider temperature range. Otherwise the Schottky 
peak of the specific heat, for instance, can not be repro- 
duced at all [19]. In this context, a mixed Bose- Fermi 
representation of spin operators [33] may be another use- 
ful scheme to thermal calculations. 

The magnetic susceptibility is calculated as 



3k n T 



{n% + 1) . 



(3.4) 



and is shown in Fig. 3, together with quantum Monte 
Carlo calculations at N = 32, which are almost the 
long-chain-limit behavior. Modified spin- wave calcula- 
tions well agree to the numerical findings, covering the 
low-temperature region to be hardly reached numerically. 
The decreasing behavior turns increasing for k^T > 3 J. 
A minimum in the susceptibility-temperature product is 
one of the most remarkable features of ferrimagnets. The 
ferromagnetic and antiferromagnetic terms n^(nf + 1) 
in Eq. (3.4), respectively, contribute increasing and de- 
creasing behaviors with increasing temperature. Since 
the ferromagnetic (antiferromagnetic) features are pre- 
dominant for S > 2s (S < 2s) [34], the present cases with 
(S, s) — (|, i) are rather biased ferromagnetically. In or- 
der to evaluate the exchange coupling constants in the 
Mn(hfac)NIT-i? systems, we further compare quantum 
Monte Carlo calculations with experimental findings. 
The susceptibility measurements on [Mn(hfac)NIT-Ph]e 
at H = 0.5 T [17] and those on [Mn(hfac)NIT-i-Pr] (X) un- 
der no field [18] are also plotted in Fig. 3, taking J/fce 
to be 370 K and 360 K, respectively. Carrying out semi- 
classical calculations, Caneschi et al. [18] estimated J/fce 
for [Mn(hfac)NIT-i?] 00 as 474.5 K (R = i-Pr), 373.3 K 
(R = Et), 311.8 K (R = Me), and 299.5 K (R = Ph). 
However, smaller values are obtained through a Fischer's 
model [35]. As for [Mn(hfac)NITPh] 6 , no quantitative 
assignment of J has been given so far. The present esti- 
mates will contribute toward establishing the standard. 



IV. NUCLEAR SPIN-LATTICE RELAXATION 

NMR measurements on isotropic magnets are nec- 
essarily performed with an applied field. An induced 
Zeeman energy gap is usually smaller than the ex- 
change interaction but larger than the nuclear energy 
scale: 10 5 Uujn ^5 10 2 </^b# ^ J- Considering the 
electronic- nuclear energy-conservation requirement, the 
Raman process should play a leading role in the nuclear 
spin-lattice relaxation [36]. The Raman relaxation rate 
is generally given by 



Ti ~ 



E, 



-E n /k B T 



-E n /k B T 



x I (m| J2( A i S i + <sl)\n)\ 2 5{E m - E n - hcj N ) , (4.1) 
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where A? and a| are the dipolar coupling constants be- 
tween the nuclear and electronic spins in the jth unit cell, 
= 7nH is the Larmor frequency of the nuclei with 7n 
being the gyromagnetic ratio, and the summation ^ n is 
taken over all the electronic eigenstates \n) with energy 
E n . Taking account of the significant difference between 
the electronic and nuclear energy scales, the relaxation 
rate (4.1) is expressed in terms of modified spin waves as 



_L 2fc(g/i B 7N) 2 



N 



<T=± 



x \W(k) - air{k)\ 2 nl{nl + l)p%k) 

k 

+ £'l^~ CT (fc) - a^ k r(k)\ 2 nUnk + l)p°(kj\ , (4.2) 



fc=0.7T> 



where ^' fe denotes the limited summation J2k — S 
ip a (k) = (S + s)/2[(S-s) 2 +4Sssm 2 (k/2)} 1 / 2 + (j/2, and 
Al and a z k are the Fourier transforms of the coupling con- 
stants, whose k dependences are hereafter assumed to be 
negligible. A contrast between zero and one dimensions 
lies in the spectral density p a (k), which originates in the 
energy-conservation requirement 8{E m — E n — Filon) in 
Eq. (4.1). In the thermodynamic limit N — > oo, p a (k) 
is definitely the differential coefficients of the dispersion 
relations, while for small clusters, it is approximately re- 
placed by the difference quotients: 



P ± (k) 
P ± (k) 



[d^(k)/dk] k=ko 
2n/N 

\uj ± (k)-uj ± (k-2Tr/N)\ 



for N -> oo , (4.3a) 
for N = 0(1) , (4.3b) 



where ko as a function of k is given by to (ko) 



= (*)■ 



ftojN = 0. If we further process Eq. (4.3a) assuming the 
predominance of k ~ contributions [37] in integrating 
Eq. (4.2), which is well justified unless temperature and 
an applied field are sufficiently high and strong, respec- 
tively, we obtain an expression 



p ± (k) 



1 



2vyjk 2 + huj-^/v 



for N — > oo. 



(4.4) 



where v = [Ss - (S + s)Ji + VSsr 2 ]J/2(S - s) is the 
curvature of the dispersion relations at small momenta. 
Equation (4.4) is in contrast with Eq. (4.3b) in that it 
depends on an applied field. 

Another consideration should be directed to the mod- 
ified spin-wave scheme in calculating the relaxation rate. 
For isotropic ferrimagnets with gapless excitations, the 
constraint (3.1) works so well as not only to suppress 
the thermal divergence of the boson number but also to 
give a precise description of the low-temperature ther- 
modynamics [38]. On the other hand, once a field is 
applied and a gap A opens in the electronic energy spec- 
trum, the boson number should exponentially decreases 
as cx e~ A / kliT at low temperatures, whereas the con- 
straint still keeps it finite even at k^T <C gp^H . Then 
we adjust Eq. (3.1) to the present situation as 



-*)££ 



= (S + s)Ne~ A/kBT . (4.5) 



This condition smoothly turns into Eq. (3.1) as H — > 
and the modification is essentially restricted to the suffi- 
ciently low-temperature region k^T < A < 10 -2 J. 
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FIG. 4. The interacting-spin-wave calculations of the temperature dependence of the nuclear spin-lattice relaxation rate with 
varying location of the probe nuclei at N — 6 (the upper three) and N — oo (the lower three). 
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FIG. 5. The interacting-spin-wave calculations of the field dependence of the nuclear spin-lattice relaxation rate with varying 
location of the probe nuclei at TV = 6 (the upper three) and N — oo (the lower three). 
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FIG. 6. The interacting-spin-wave calculations of the nu- 
clear spin-lattice relaxation rate as a function of inverse square 
root of an applied field with varying location of the probe nu- 
clei at N = 6 (the upper two) and N — oo (the lower two). 

Now the relaxation rate is calculated for (MnNIT)6 
and (MnNIT)oo. Besides temperature and an applied 
field, it may be a function of the location of probe nuclei, 
which can be described by a parameter r = A z /a z ~ 
(d s /ds) 3 , where ds and d s are the distances between the 
nuclear and electronic spins. Figure 4 shows 1/Ti as a 
function of temperature. The ferromagnetic and anti- 
ferromagnetic spin waves contribute different tempera- 
ture dependences to 1/Ti and their weights vary with 
the parameter r. Considering the predominance of the 
k ~ contributions in Eq. (4.2), at r ~ ■0 <T (O)/-0"' T (O) = 
(S/s)' 7 , the a excitation mode hardly mediates the nuclear 
spin relaxation, in other words, the electronic excitations 
of a mode are invisible to the nuclear spins. The temper- 



ature dependence of 1 jT\ is indeed of antiferromagnetic 
(ferromagnetic) aspect at r = 1/5 (r = 5), while other- 
wise it is of mixed aspect. If the probe nuclei are located 
as r ~ 1/5, the relaxation rate is extremely small. These 
calculations can be observed by taking different kinds of 
nuclei as probes, such as 1 H, 13 C, 19 F, and 55 Mn in the 
present systems. Equation (4.2) is still valid for r — > oo, 
which corresponds to 55 Mn NMR. 

Figure 5 shows 1/Ti as a function of an applied field. 
There appears a clear contrast between zero and one di- 
mensions, where 1/Ti 1S saturated and diverging, respec- 
tively, with decreasing temperature and field, as long as 
the ferromagnetic spin waves are visible to the nuclear 
spins. In finite clusters, any field dependence of 1/Ti 
is necessarily attributed to n k , which are simply expo- 
nential with respect to H, whereas in long chains, 1/Ti is 
more varied with a field, depending on it through both h k 
and p ± (k). Therefore, unless the Zeeman energy becomes 
comparable to the exchange interaction, 1/Ti exhibits lit- 
tle field dependence in finite clusters. In order to bring 
out their characteristic field dependences more quantita- 
tively, we plot 1/Ti as a function of 1/yH in Fig. 6. 
nZ and p~{k) are both peaked at k = 0. More and more 
weight centers on k = with decreasing temperature and 
field for N — > oo in particular. Therefore, at low temper- 
atures, the N — > oo k integration in Eq. (4.2) may ap- 
proximately be replaced by the k — contribution, which 
is in proportion to 1/s/H. Thus, as long as an applied 
field is moderate, a 1/yH -linear behavior is observed at 
low temperatures. With increasing temperature it turns 
logarithmic, sloping more gently, due to the /c-integration 
effect. Under strong fields, they are all masked behind the 
overwhelming Zeeman term oc e~9^BH/k B T com i n g f rom 
n^T . At r = 1/5, where only the antiferromagnetic spin 
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waves are active for 1/Ti, the nuclear spins exhibit ex- 
tremely slow dynamics. Since is not peaked even at 
low temperatures, there appears no l/\/H dependence. 
With increasing temperature and field, the antifcrromag- 
netic excitation branch lowers in energy and thermally 
assists the relaxation, ending up with increasing 1/Ti. 
A minimum ofl/Ti as a function of H can be observed 
only when the ferromagnetic spin waves are almost off. 
In finite clusters, the discrete spectrum may in principle 
lead to 1/Ti oscillating as a function of H. 



V. CONCLUDING REMARKS 



Motivated by inorganic-organic hybrid compounds, 
[Mn(hfac)NIT-Ph] 6 and [Mn(hfac)NIT-i?] oc , we have 
demonstrated model calculations of the low-energy spin 
dynamics in ferrimagnctic clusters and chains. Tempera- 
ture dependence of 1 /T\ drastically varies with the loca- 
tion of the nuclei in both compounds, though the relax- 
ation time scale is much larger in zero dimension than in 
one dimension. There are special points for the nuclei, 
characterized as (ds/dg) 3, ~ s/S, where the nuclear spin 
relaxation can hardly be assisted by the low-lying ferro- 
magnetic excitations of the electronic spins and therefore 
extremely slow dynamics is observed. 

Field dependence of 1/Ti in long chains forms a strik- 
ing contrast to that in finite clusters, diverging with 
decreasing field at low temperatures. The present ob- 
servations should be distinguished from the 1/VH or 
\n(l/H) dependence of diffusion-dominated dynamics 
[39,40] , which originates from transverse spin fluctuations 
and distinctly appears at high temperatures. In the vicin- 
ity of the special points of (d s /ds) 3 ~ s/S, a minimum 
of 1/Ti as a function of H can be observed. 

Besides the Mn(hfac)NIT-i? systems, there are 
a series of ferrimagnctic bimetallic chain com- 
pounds MCu(pba)(H 2 0) 3 -nH 2 (M = Mn,Ni; 
pba = 1, 3-propylenebis(oxamato)) [41] and 
MCu(pbaOH)(H 2 0) 3 -nH 2 (M = Fe,Co,Ni; pbaOH = 
2-hydroxy-l,3-propylenebis(oxamato)) [42], which are 
also describable by the Hamiltonian (2.1). Since their 
exchange coupling constants are much smaller than those 
of the Mn(hfac)NIT-i? systems, they are complemen- 
tary in field-applied measurements. We hope that the 
present calculations will stimulate further experimental 
explorations into quantum dynamics on the way from 
zero- to one-dimensional magnets. Nuclear spin-lattice 
relaxation-time measurements on [Mn(hfac)NIT-Ph]g 
and [Mn(hfac)NIT-i?]oo are strongly encouraged. 

This work was supported by the Ministry of Educa- 
tion, Culture, Sports, Science, and Technology of Japan, 
and the Nissan Science Foundation. 
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